The cosmological background of vector modes 
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We investigate the spectrum of vector modes today which is generated at second order by density 
perturbations. The vector mode background that is generated by structure formation is small but 
in principle it contributes to the integrated Sachs- Wolfe effect, to redshift-space distortions and to 
weak lensing. We recover, clarify and extend previous results, and explain carefully why no vorticity 
is generated in the fluid at second order. The amplitude of the induced vector mode in the metric 
is around 1% that of the first-order scalars on small scales. We also calculate the power spectrum 
and the energy density of the vector part of the shear at second order. 



I. INTRODUCTION 

Cosmological vector perturbations at linear order satisfy an evolution equation and a momentum constraint equation 
which can be found by calculating the i — j and — z parts of the Einstein field equations respectively. In the case of a 
perfect fluid (or scalar field) there is no source in the vector evolution equation and its solutions decay as 1/a^. In the 
standard models of structure formation, inflation does not generate vector perturbations, and we can therefore ignore 
vector perturbations at flrst order. However, vector perturbations will be generated at nonlinear order by the growth of 
density perturbations. Here we revisit the analysis of this cosmological vector background (see Refs. [H |2 [2 HI El E] 
for previous work). As we show below, the amplitude of the induced vectors is surprisingly large on small scales, when 
compared to the amplitude of the first-order Newtonian potential. In order to investigate whether this is important 
dynamically we calculate the vector part of the shear and compare it to the usual first-order scalar shear. 

We consider perturbations of a flat Robertson- Walker background up to second order: g^^, — gfj,v + S'-^^^ gfj,u + S^^^ g^j^^. 
At first order, S'^^^gf^i, only contains scalar perturbations: we neglect the tensor perturbations, and vector perturbations 
are not generated at first order in the standard model. For the second-order perturbations, S^'^^g^^,, we project out 
the vector modes. In the Poisson gauge [T] 



[1-H2$(i)] drf 



a^S^f''dx'dri- 



[l - 2$(i)] dx' 



(1) 



where '^(i) is the first-order Newtonian potential (we assume zero scalar anisotropic stress at first order), and S^^^) 
describes the gauge-invariant second-order vector modes, so that cJiS'^g) = 0. We effectively ignore the second-order 
scalar mode <I>(2) since we are interested only in the second-order vector modes: second-order modes can be consistently 
split into scalar and vector, but the equations for the second-order vector modes will contain source terms that are 
quadratic in the first-order scalar modes. In what follows we will drop the order indices when there is no ambiguity. 
With this in mind, the fiuid four- velocity is given by 
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wher 



e wp'' = diV(i) and 9iW(2) — 0. 



The background dynamics are given by 



-a^A, n' = -ATiG{l + w)a^p + T-l^ 



(3) 
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where w = p/ p = const (later we specialize to w = 0). The first-order perturbed field equations lead to [5] 

"('^ " 47rGpa2(l + u.) ' ^ ^ 

' 

where the Newtonian potential obeys the evolution equation 

+ m{l + + [(1 + w)W - wV^] $ = . (6) 

All linear scalar modes are determined by <i>, and terms quadratic in $ and its derivatives will source the second- 
order vector modes. Therefore we will require the power spectrum V<^. In the matter era (w = 0), neglecting A, we 
have a = 00(77/770)^, and the solution in Fourier space is <i>m(fc, 77) — Amik), where we remove the decaying mode. 
Then the power spectrum today is computed using the growth suppression and transfer functions. (See Appendix A 
for more details.) 

Vector perturbations typically produce vorticity and a transverse shear in the fluid four- velocity, Eq. ^ . In order 
to compute them, we need the covariant definitions (see [S] for a recent review): 

w^j. = S^^'^ujfj^^ + S^'^'^ujfj^.y = hf^h'^^Ua-iB , (7) 

= S^^^a^, + S^^^a^, = l^hfX) - "a;/? , (8) 

where h^i, = g^v + u^u^, is the projector into the instantaneous fiuid rest space, and ui^iuu'^ = = a^^u" . The vorticity 
is a purely vector quantity. The vector part of the second order shear is defined via (see Appendix D for how to get 
Gj from (^^^Vij ) 

(5(2)^,, = a%a,) , 9,(7^ = 0. (9) 

At first order, there is only scalar shear: 

5^^^uj,j = 0, (10) 
<5(iV,, = a (d,dj ~ ^^y V^") v(^) . (11) 



II. SECOND ORDER VECTOR MODES 



A. Vorticity 

We use a covariant and fully nonlinear approach in this sub-section, which is more direct and transparent than a 
perturbative approach in the case of vorticity, and which also leads to a more general result. 

The vorticity tensor of a fluid with four- velocity deflnes the vorticity vector — e^^^uo^'^ /2, where e^,^-y is the 
covariant permutation tensor in the fluid rest space. The vorticity vector obeys the following propagation equation, 
which is covariant and allows for full nonlinearity [9; : 
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h^J'Lo^ ^ --Qijj^j^ - -cw:\u^ + a^i,Lo'' . (12) 

Here Q — V^w'^ is the volume expansion rate, an overdot denotes covariant differentiation along the fluid flow (u^V^), 
so that iifj, is the fluid four-acceleration, and the covariant spatial curl is defined by curln^ ~ e^i^^D'^n''', where is 
the spatially projected covariant derivative (D^rij/ = hp^^hiy^V anp) . 

In order to evaluate the curl of the acceleration, we need the momentum conservation equation [5], 

4 

hf.'^q^ + -6g^ + {p + p)ut^ + ^^LP + D^TTpj. + CTt^uQ" - e^^jUJ^q'^ - TTfj^^ii" = , (13) 



where is the spatial momentum density fiux relative to and tt^^ is the fiuid anisotropic stress (spatial, tracefree 
and symmetric). In general, Eq. (13 1 shows that the curl of acceleration can introduce source terms for vorticity in 
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Eq. (12 1. In the case of a perfect fluid however, we have 
sound speed. Then using the exact identity [9] 



= nfj_^ and D^p 



curlD„/ = -2/w. 



c^D^p, where Cg is the adiabatic 



(14) 



the curl of Eq. (13) gives (p + p)curl'u^, = 2p uj + e ^^-yii'^'D'^ (p + p) . Using Eq. (13l again, we find that the second term 



is proportional to e^^-^D^pD^p , which vanishes, since D^p is parallel to D^p for a perfect fluid. Collecting results, we 
arrive at the fully nonlinear vorticity propagation equation for a perfect fluid: 



(15) 



(Note that the last term on the left is at least third order.) This equation shows that there is no source for vorticity, 
so that vorticity cannot be generated in a perfect fluid, at any perturbative order. 

In particular, there is no generation of vorticity at non-linear order by first-order scalar perturbations, in the case 
of a perfect fluid, and thus 



0. 



(16) 



If there is primordial vorticity, then it must be introduced as an initial condition. Any primordial vorticity will 
simply redshift away as the universe expands, according to Eq. (151, and will be entirely unaffected by the growth 
of density perturbations. Effectively, the density perturbations generate metric vector perturbations, and the fluid 
velocity adjusts so as to maintain zero vorticity. This is similar to what happens with the Harrison mechanism for 
magnetogenesis, where the vector modes generated at first order by defects cannot induce vorticity in the plasma |10j . 
For more than one perfect fluid, vorticity non-generation applies separately to each fluid, as long as there is no 
momentum exchange between the fluids. This is the case for example with cold dark matter and baryons, which 
interact only gravitationally. 



Returning to the perturbative analysis, we can now use Eq. ( 16 1, together with the Einstein equations, to determine 
the vector metric perturbation S"^. 



B. Vector metric perturbations 

The second-order vorticity is given by Eqs. ([2]) and (0: 



(17) 



which is in agreement with the expression in [11 . Then Eq. ( |16[ ) becomes a constraint on t)^2) ~ •S'*- However, we have 
another constraint from the second order Oi Einstein equation [7] 



,(2) 



s, 



= -V^S*, + 8 



2$'9»$ + ^^^2j^^'* 9, ($' + H$) 



(18) 



where V denotes schematically the vector part of the quadratic source term which can be extracted in Fourier space 
as shown in Equation (C2|. Substituting Eqs. (|4| and (18 1 into Eq. (17 1, we obtain 



Then Eqs. I\16] and ([19f imply that 
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— — |v2$a,($' + 7i$)| 

1 + ui) L J 



(19) 



(20) 



This directly recovers the solution that was obtained in Ref. [7] via a Fourier-space projection of the second-order ij 
Einstein equation, i.e., the vector anisotropic stress constraint. 



(21) 
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Here ttj (with djir^ = 0) is the effective vector anisotropic stress from second-order density perturbations (see [7] and 
Appendix B). 

To evaluate Si via Eq. (20 1, we work in Fourier space, and the details are given in Appendix|C] Using the first-order 
power spectrum (Appendix A), the resulting second-order vector power spectrum Vs is shown in Fig. [l] This may 
be compared with the power spectrum in [7 , which is computed from scalar modes in the radiation era (w — 1/3), 
assuming a power-law form for the scalar spectrum. Here we have found Vs today, and we have calculated 7^$ directly 
from the first-order solutions, using the transfer function to relate back to the primordial perturbations and the growth 
suppression factor to take account of A in the background. 
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FIG. 1: The power spectrum today of metric vector modes generated at second order by density perturbations. Left: Vs 
together with the power spectra of first-order quantities: the density perturbation 5 = Sp/p and comoving density perturbation 
A = S-3Hv'-^\ and $ (using best-fit WMAP5 parameters [T^]: fi^/i^ = 0.1326, lOOfJt/i^ = 2.263, /i = 0.719). Right: We show 
how increasing the baryon fraction decreases the power in the vectors above keq. The baryon oscillations are washed out to 
some extent in the vectors as can be seen by comparing Vs with the first-order transfer function in the top panel. 



In Fig. [Tjwe compare the power spectrum of the vectors with that of the Newtonian potential. For comparison we 
also show the power spectrum of the two gauge-invariant density perturbations. The amplitude of the vectors decays 
on small scales, k > k^q ~ 0.009 Mpc~^, in contrast to the density perturbation, which is growing, but in line with <f>. 
For WMAP5 data [E] 

Vs ~ 6.5 X IQ-^V^ for k > /c^uk ~ 0.09 Mpc"^ , (22) 

so that the amplitude of the metric vector modes is nearly 1% that of the metric scalar modes on small scales. For 
zero baryons we find that Vs ~ 2;-i(5.49f}„/i2 - 0.13f -^^V$ ~ {Inkf/k^ for k > fcsuk « 0.09 Mpc"^ On large 
scales Vs scales like k, with a peak in the spectrum around the equality scale. This is analogous to the peak in the 
induced gravitational wave background on similar scales |13j . In Fig. [Tjwe also show that including baryons induces 
oscillations in the vector power spectrum (bottom panel). They are washed out in comparison to the those present 
in the scalars (top panel), but are still very prominent. 

The overall shape of the S*' spectrum may be understood from the generation of vectors during the radiation 
era [7]. Vector modes grow outside the Hubble radius as a^/^ only through the interaction of scalar modes which are 
larger than the Hubble radius. Inside the Hubble radius, vector modes decay, shghtly less rapidly than a^^, when 
fluctuations in the radiation fluid no longer support vectors. At the end of the radiation era, vector modes with 
k < keq have acquired a tilt because modes are more aggressively produced by scalars which are close to the Hubble 
radius - very long wavelength modes interact only weakly. After equality, all vector modes grow at the same rate, so 
that those which entered the Hubble radius before equality are suppressed. 
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Vector Shear 



The second order vector shear follows from Eq. (|8| as 
Using Eqs. (|3]), (|4]) and ([18|, we have 



(23) 



2a 



m^n^ (1 + «;)^{^' + " ""^^^ + ""^^ + ^^^^ + iw)^^^'^,^' 

2re [1 + 3(1 + an)(l + (f^^Sj*)' - (1 + 2u;)V^%$9j)($' + H$)}^. 



(24) 



In order to find the vector power spectrum V^v for Gi, we substitute Eq. (20 1 into Eq. (24 1, and then apply the 



vector extraction operator in Fourier space. Further details are given in Appendix C. The resulting power spectrum 
has a similar shape to Vs, since ai = Si/2 + small corrections by Eq. (24 1. 



It is also useful to compare the first-order scalar and second-order vector contributions to the shear energy density. 
We define the dimensionless shear density 



(25) 



The scalar and vector contributions to shear are given by tr = 'y(i) [Eq. (Ill] and di [Eq. (24)] respectively. They 
define scalar power VaS and vector power Va-v, which then define the spatial averages of fl^s^ ^a-v via 



dink 



dink 



12^2 



(26) 



Note that there are 2 polarizations implicit in VaV- 

The quantities in Eq. (26l are shown in Fig. [2] As k 0, d{fla-y) /dink ~ /c'^ and d{VL„s) / dink ^ k*. This shows 
that the second-order vector shear is much smaller than the linear scalar shear - except on very large scales. However, 
this feature in the Poisson gauge will not lead to any growing physical effect, since both quantities are decaying on 
large scales. 
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FIG. 2: The energy densities of the first-order scalar and second-order vector shear per logarithmic fc-interval. 



The vector part of the shear may be interpreted as a rotational quantity, even though the vorticity of the fluid is 
zero. We can see this readily via the covariant approach of section IIA. The vector shear is 



0. 



(27) 



6 



By boosting from the fluid rest frame to a different frame, u'^ = u'^ + (with v^^u^ = 0), the vorticity and shear 
become [9j 

= - ^curl Vf^, = (Tf^ + . (28) 

If we choose the frame by = — c^, the vector part of the shear is zero, but the vorticity o)^ no longer vanishes. 
This shows the essentially rotational nature of the fluid vector shear, even though the vorticity of the fluid is zero. 
Note that it is not possible to boost away the scalar or tensor part of the shear in this way. 



III. CONCLUSIONS 



We have computed the power spectra for the metric vector perturbations and vector shear at second order, generated 
by first-order scalar perturbations in a ACDM model. In addition, we used a covariant approach to show explicitly 
how vorticity is not generated in a perfect fluid at any perturbative order by first-order scalar perturbations, so that 
there is no vorticity in the matter. In order to generate vorticity, one requires either an imperfect fluid, or momentum 
exchange between the fluid and another fluid. Momentum exchange (via Compton and Coulomb interactions) between 
electrons, protons and photons in the radiation and recombination eras can generate vorticity and magnetic flelds at 
second order HH [HITil [H [T71 [HI . 

In order to obtain the vector quantities and vj^^) 
constraint [Eq. (18l]. Alternatively, one could also use the ij 
conservation equation, which has the form 



we used the vanishing of vorticity [Eq. ( 16 1] and the Oi Einstein 
Einstein equation [Eq. (21 1] and the momentum 



{vi2)^ - + H{1- 3cl) {«(2), -S,}^M,, 



(29) 



where the source term is given in Appendix B. 

The cosmological background of vector modes is small, especially if measured in terms of the dimensionless shear 
density, as shown in Figs, [l] and |2] However, given that the amplitude of the vector modes in the metric is as large as 
~ 1% of the metric first-order scalar modes, in principle these vector modes will have an effect on various cosmological 
observations. In particular: 

• Redshift-space distortions |20]: the divergenceless velocity w^j) "^i^l make a contribution to radial peculiar 
velocities and thus to redshift-space distortions. 

• Large-angle CMB temperature anisotropics \T, : the vector modes will contribute to the Doppler and integrated 
Sachs- Wolfe effects: 



1 



(30) 



• Weak Icnsing [2ll |22] : vector modes produce a deflection angle 

dx(v X s) X 

E ^ ' 



(31) 



• CMB polarization [31[53]: the vector modes will leave a characteristic imprint on CMB polarization. 

Further work is needed to compute the size of these vector corrections. They arc likely to be signiflcant mainly 
below ^ lOMpc, but in this regime the scalar non-linear effects are important and likely to dominate. 

The vector degree of freedom forms an integral part of the perturbative expansion when one goes beyond linear 
order. At the order we have considered, vectors must be present essentially through a constraint in the field equations 
arising at order <i>^, even though vectors have no independent propagating degrees of freedom. This is distinct from 
the intrinsically propagating degree of freedom in the scalar-induced gravitational wave background [13i i24i i25j . We 
have shown that the vector mode background has maximum power around the equality scale, similar to the induced 
gravitational wave background, and the metric vector modes achieve their maximal fraction of the linear metric scalar 
modes on scales below the Silk scale. The size of the vector contribution to the full non- linear power spectrum relevant 
for structure formation remains to be calculated. 
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APPENDIX A: FIRST-ORDER SCALAR POWER SPECTRUM 



The power spectrum for the first-order scalar perturbations is defined by 

(<i>*(fe, ,7)<i>(fc', r^)) = ^<53(fc - k')V^{k, 77). 



In the early radiation era. 



^3 91 Rty"^ 



4867r2 



(Al) 



(A2) 



where is the primordial power of the curvature perturbation, with [12] « 2.41 x 10~^ at a scale kcMB — 
0.002Mpc~^. By conservation of the curvature perturbation, = 9<i>r/10 at equality, so the early matter power is 
given by 



Am{k) — — Ar{k) « — 



(A3) 



assuming a scale-invariant initial spectrum. Here ijoo is a normalization parameter in the A growth suppression 
function 



1 + ^^m{z) 



(A4) 



and goo is chosen so that g(0) = 1. The power today is given by A„(fc) = A„i(k)T{k), where T(fc) is the normalized 
transfer function [26j (see Fig. [T|): 



/3Ak 



g'T{k)' 



(A5) 



APPENDIX B: VECTOR EVOLUTION AND MOMENTUM CONSERVATION 



The source term in the vector evolution equation Eq. (21 1 is as follows (where we assume w = const) 



a9(i7rj) = 4 



^^{9,$'5,<I>' + 2H%$9,)$'} , (Bl) 



where V denotes the vector part, and we have used the fact that {<i>d,dj<i>} = -{^^^^j<^>} . 

In Eq. (29), the source term in the vector momentum conservation equation is as follows (where we also assume 
w — const): 
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9H^f2^(l + wfM, = 6n^ [24(1 + clf - 12(1 + cl){2 + w) + n„^{l + w){5 + ISc^ - 15w)] {$9,$'}^ 
+6H^ [24(1 + clf - 12(1 + cj){2 + w) + 2n„,{l + w)(3c^ - 2)] {$'9,$}^ 

+2n [10(1 + cf) - 18(1 + clf - 6(1 + c2)(i + - 4] {Ha,$v2$ + a,$'v2$}^ 
-36^2 c2 [2(1 + c^) + n„,{i + w)] {^v^a,*}'' + 8cl {Hv2$'a$ + v2$'ai$'}^ 
+4(1 + cl)cl {v^^v^a^*}^ - 247i:(i + cl)cl {$'v^a$}^ 

+8H {d'^^'dkd^^ + d''<i>dkdi^'}^ . (B2) 



APPENDIX C: VECTOR MODE POWER SPECTRUM 

We define the Fourier transform of the vector perturbation as 

S,{x, r?) = ^^-i^ J d^k [S{k, ri)e,{k) + S{k, ,7)e,(fc)] e^^^, (CI) 

where the two orthonormal basis vectors e and e are orthogonal to k with 

S{k, ri) = J d^xS,{x, ,7)e»(fc)e-*-. (C2) 

The inverse integral is similarly defined for its second parity S(k,ri). Equation (20 1 gives a solution for S{k,ri) in 
Fourier space 

^(fc.^) = 3^^21)3/2 ^ / d'k' \k-kfk',B{k-k',k',r,), (C3) 

where 

S(fei, fe2, V) = W"'*(fcl, 77) [$'(fc2, ??) + W*(fe2, • 

Defining the power spectrum as 

27r2 

{S* {k, r,) S (fc', r,)) = (fc _ fc,) (C4) 

we find, using Wick's theorem, 

^^^^'^^^=9^^ / d'k'\k-kf[k'^e'{k)] [K^e^{k)]B{\k-klk',rj) 



x{\k- k'\^B{\k - fe'l, fc', r;) - {k'fB{k', \k ~ fe'|, ry)} 

This may be simplified to give 



where 



(C5) 



), (C6) 



nOC pV-\-l 

n(0= / dv / du ^ {uvy^iu^ ^v^)[4v'^ ^ {l + v^ -uY][T{kv)T{ku)]\ (C7) 
and V = fc'/fc, M = ^/l + v'^ — 2v cos and cosf? = fc' • k/{k'k). 



APPENDIX D: VECTOR SHEAR POWER SPECTRUM 

To extract the divergenceless vector (Ti from 9(icrj) , we use the operator [7] 



2i 



(27r) 



(Dl) 



Then (Jm{x) ~ V',y{{x,x')d(iaj){x'), and Eqs. (20l and (24) lead to the second-order vector shear in Fourier space 

-2ai e^{k) 



C{k',k~k',r])+ 1 



B{k',k~k',T]) 



where 



C(fcl, fc2, V) - ^7) ['i''(fc2, 77) + nHk2, V)] ■ 



The power spectrum is defined as in Eq. (|C4|. By Eq. (|D2| and Wick's theorem, we obtain 
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Vav{k,vi) 
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I -2.g2n(2 + u2 + 3i;2) 
n(l + 2t;2) + 



1 - ^ I 



(D2) 



(D3) 
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where g is given by Eq. (A4| and 11 by Eq. (C7). 
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